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Abstract 

The bivariate high-temperature expansion of the spin-spin correlation-function of the three- 
dimensional classical XY (planar rotator) model, with spatially-anisotropic nearest-neighbor cou- 
plings, is extended from the 10th through the 21st order. The computation is carried out for the 
simple-cubic lattice, in the absence of magnetic field, in the case in which the coupling strength 
along the z-axis of the lattice is different from those along the x- and the y-axes. It is then possible 
to determine accurately the critical temperature as function of the parameter R which characterizes 
the coupling anisotropy and to check numerically the universality, with respect to R, of the critical 
exponents of the three-dimensional anisotropic system. The analysis of our data also shows that 
the main predictions of the generalized scaling theory for the crossover from the three-dimensional 
to the two-dimensional critical behavior are compatible with the series extrapolations. 
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I. INTRODUCTION 



The three-dimensional(3D) layered magnetic spin systems in which the strength of the inter- 
actions among the layers is much smaller than within the layers, are often referred to as quasi- 
two-dimensional. Although strictly two-dimensional(2D) magnetic systems do not exist in nature, 
their statistical mechanics can be studied by experimenting with diverse more mundane structures 
and in particular by exploring how quasi-two-dimensional systems^ crossover, i.e. change their uni- 
versality class, on going from the 3D to the 2D critical regime. From a more general standpoint, 
the study of spatially anisotropic systems also provides the simplest example of a wide variety of 
crossover phenomena£>^£>£ of different origin which may occur near criticality. 

The simplest Hamiltonian which can model a quasi-two-dimensional magnetic system, in the 
absence of a magnetic field, is that of an iV-vector spin model with axially anisotropic couplings 

H an {v} = -iVJi Y, $(f) ■ ~ NJ 2 *(f) ■ v(r ') (1) 

nn(xy) nn(z) 

We have indicated by v(r) an TV-component classical spin-vector of unit length located at the site 
f of a simple-cubic (sc) lattice. The first sum in ([1]) is extended to the nearest-neighbor (nn) spin 
pairs within each horizontal (xy) layer, while the second sum is over the nn spins in adjacent layers 
along the ^-direction. We shall denote by R = J%j J\ the ratio of the interlayer to the intralayer 
coupling strength which characterizes the spatial anisotropy of the spin couplings and therefore is 
sometimes referred to as anisotropy parameter. The thermodynamic quantities of the model can 
be expressed as functions of the variables K\ = J\jkT and K2 = J2/kT, with k the Boltzmann 
constant. One may, equivalently, choose either the pair of variables K\ and R or the pair K2 and 
R = 1/R. For R — > 0, the system becomes a stack of non-interacting spin layers. When R = 1, the 
system has directionally isotropic interactions. For R — > 00, or equivalently for R — > 0, it reduces 
to an array of non-interacting spin chains. 

Only a few pioneering studies^ ^ 10 ' 11 ' 12 of the Hamiltonian ([1]) by high-temperature(HT) meth- 
ods are presently available. They were aimed at: 

1) a numerical test of the critical universality^ for the anisotropic system, in particular of the 
i?-independence of the critical exponents as long as R > 0; 

2) an investigation of the change of universality class of the critical transition as R — > 0, i.e. of 
the crossover from the sc to the square-lattice critical behavior. 

They relied on HT series expansions, in terms of the two variables K\ and K2, computed 
through 11th order 8 for N = 1 (the spin 1/2 Ising model) and through 10th order— for N = 2 
(the planar rotator or XY model) or N = 3 (the classical Heisenberg model), on the sc lattice. 
The corresponding expansions for the face-centered-cubic lattice also reached the same orders. 
Altogether 78 coefficients were computed in the sc-lattice Ising case and 66 coefficients in the other 
cases but, unfortunately, no higher order coefficients were added since. In the Ising case, although 
rather short, the expansions are sufficiently well behaved that their extrapolations can support 
unambiguously the simplest theoretical expectations concerning the crossover. The accuracy of the 
first series analyses in the Ising case could be only marginally improved— by resorting to bivariate 
Partial-Differential approximants^ii^ (i.e. by approximately resumming the HT expansions in 
terms of the solution of a linear first-order partial differential equation with appropriately chosen 
bivariate polynomial coefficients) instead of using the conventional ratio or Pade approximant 
(PA) methods. The reason of this failure is that there is no substitute for significantly longer 
expansions. Later on, also several MonteCarlo simulations^ 1 ^ were carried out in an attempt at 
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further clarifying the crossover behavior of Ising systems, but the accuracy of the results is still 
subject to controversy for small positive R, i.e. in the region of main interest. On the other 
hand, for the models with N > 1, the earliest HT analyses were inconclusive, thus calling for a 
substantial extension of the expansions. In particular, the crossover issue remained to be studied, 
because even the existence of a critical point at nonzero temperature and the nature of the critical 
singularity in the 2D limit were not yet firmly assessed at the time of the first analyses. Also the 
successive simulation studies, carried out when the critical behaviors in 2D of the N = 2 model^ 
(and of the N > 2 models 2 ^) were better understood, probably cannot yet be considered sufficiently 
accurat e 21 1 22 or are not directly comparable 2 ^ with the series analyses. 

We have been motivated by this situation to devote our study to the N = 2 anisotropic model, 
taking advantage of our recent extensions through order 21, (i.e. from 66 to 253 series coefficients), 
of the bivariate HT expansions for the two-spin correlation-function and its moments on the sc 
lattice. Another reason of interest into this model is that it has been suggested to provide an 
approximate description of layered high-T c superconductors^. 

The paper is organized as follows. In Sec. II, we briefly mention the algorithm adopted and 
specify the results of our series computations. In Sec. Ill, we discuss some of the simplest predictions 
of the extended phenomenological scaling theory for the crossover from the 3D to the 2D (and 
from the 3D to the ID) critical behavior. We begin by reviewing in some detail the well studied 
N = 1 case only to recall the general ideas of this approach and to contrast its features with those 
of the less studied and more complex N > 2 cases. In Sec. IV, we outline our numerical analysis of 
the expansions and compare the results with the theoretical expectations. 

We should finally notice that, throughout this report, for reasons of clarity, we have used a 
notation sometimes different and heavier, but more detailed and perhaps more explicit, than that 
generally adopted in the earliest studies of the crossover phenomena. 

II. EXTENDED HIGH-TEMPERATURE EXPANSIONS 

Our computation of the series coefficients was carried out using a computerized recursive al- 
gorithm based on the Schwinger-Dyson equations 25 . This method was initially applied only to 
the determination of single-variable HT expansions. Only recently, taking advantage of the great 
improvements of the computer performances of the last decades, it could be straightforwardly 
adapte d 26 ' 27 to derive also the more memory-demanding and computationally-intensive bivariate 
expansions for a wide class of isotropic and anisotropic XY models with nn and next-to-nn in- 
teractions. It should be noted that in our approach only extended-integer exact arithmetic is 
used, thus avoiding all roundoff errors which limited the precision of the preceding^ floating-point 
computations. To give an idea of the performance of the algorithm, let us note that an ordinary 
single-processor desktop personal-computer(pc) can complete, in less than a second, all lOth-order 
calculations for the anisotropic N = 2 case so far documented^ in the literature. The calculation 
of the next nine orders takes a few days. To compute the last two orders, we have used a pc-cluster 
for a time equivalent to approximately six months of a single pc. 

Fixing N = 2, we have calculated the spin-spin correlations 



for all values of x for which the HT expansion coefficients are non-trivial within the maximum 
order reached. As usual, here < O >= Tr(Oexp[—H an ])/Tr(exp[—H an ]). 




(2) 
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In terms of (|2|), we have formed the expansions of the Zth order spherical moments of the 
correlation- function: 

m®(N;K 1} R) = \x\ l < v(0) ■ v(x) > (3) 

X 

and, in particular, of the reduced ferromagnetic susceptibility defined as x(N; Ki> R) = 
m^\N;K u R). 

The second-moment correlation- length is expressed, in terms of m^ 2 \N; Ki, R) and x(iV; K±, R) 

as 

[((N-^Rjf = m^(N;K 1 ,R)/2d X (N;K 1 ,R). (4) 

where d is the lattice dimensionality i.e. d = 1 for R = 0, d = 2 for R = and d = 3 otherwise. 

Usually, for the univariate HT expansions, the only accessible validation procedure of an ex- 
tended computation is the comparison with the lower-order results that might be already known. 
In the case at hand, only the expansion coefficients of the second moment of the correlation- 
function are tabulated through 10th order in Ref . [IcI] . but they contain small roundoff errors in 
the 8th figure at highest order. After correction of these errors they agree with our results. Of 
course, this is not a very stringent test of correctness. However, the extended expansions can be 
subjected, at all orders, to additional tests, some of which deriving from equations of Section III. 
First, we can check that, taking R = 1, the coefficients of the single- variable expansions for the 
corresponding quantities of the isotropic 3D XY model, already known2£ through order 21, are 
reproduced. Moreover, we can observe that for R = and R = 0, the expansions of x{N;Ki,R) 
and m( 2 )(iV; Ki,R) reduce, as they should, to those of the corresponding quantities in the 2. D 29 i 30 
and the ID XY model, respectively. Finally, we take advantage of eqs. (fT5|) and (f23j) of the next 
Section, in the case of the susceptibility, (or eqs. (|22p and (|24p for the second moment), to pin down 
two more among the r + 1 series coefficients occurring at the rth order. The success of this variety 
of tests, through all orders we have computed, strengthens the confidence that our extensions of 
the bivariate expansions are correct. 

Our series data for the nn correlations, the susceptibility and the second moment of the corre- 
lation function are tabulated in an appendix, which for editorial reasons was not included in the 
printed version^ of this paper. 

III. DIMENSIONAL CROSSOVER 

For all values of N, at fixed R > 0, the 3D spin models described by the Hamiltonian eq.([T]) 
display a conventional power-law critical transition. As the reduced deviation t(N; R) = 1 — 
K\/K\ C (N]R) from the critical point K\ C (N\R) of the 3-D system with anisotropy parameter R 
tends to zero from above, for the susceptibility one has x(N;Ki,R) ~ [t(N; R)]~^ N ' R \ while 
for the correlation-length one has £(N;Ki,R) ~ [t(N; R)]~ U ^ N ' ,R \ The universality hypothesis 13 
dictates that, for a given value of N, the critical exponents 7(iV; R) of the susceptibility and 
v(N; R) of the correlation-length of the AT-vector system with arbitrary finite anisotropy should be 
independent of R as long as R > and thus should coincide with the exponents 7(-/V; 1) and v{N\ 1) 
of the isotropic system. The initial part of our analysis of the HT expansions in the subsections A 
and B of Section IV, will be devoted to the determination of the critical temperature and exponents 
as functions of R for R > 0, thus making it possible to test numerically the universality of the 
exponents with respect to R. 
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When R — » the sc lattice system crosses over to a stack of uncoupled square-lattice systems and 
we expect that an anomalous behavior at criticality indicates a discontinuous change of universality 
class. The crossover behavior is described by a phenomenological scaling theory introduced in 
Ref.fl and subsequently extended and clarified in Refs.jjjajg]. This approach, whose validity has 
been verified in the mean-field approximation and in the spherical model^, is first outlined for 
N = 1 (the Ising model) in the subsection A and then generalized to cover also the case N > 2, in 
the subsection B. The predictions of the scaling theory for N = 2 will finally be compared to the 
HT-based approximations in the subsection C of Section IV. 

An anomalous behavior is also expected to occur in the R — > limit in which the sc lattice 
crosses over to the linear lattice. 

A. The N = 1 model 

For N = 1, both the 3D and the 2D TV- vector models display a power-law critical behavior. As 
a consequence, for all critical exponents 7(1; R), v(l; R), . . ., of the 3D Ising model with anisotropy 
R, the R — > limit exists and yields the corresponding exponents 7(1; 0), u(l;0),... of the 2D 
Ising model. Thus, for R = 0, one can write 

x(i;Ki,0)«Xas(i;^i,0)~r-^°) (5) 

and 

e(l;^i,0) «U(l;^i,0)~r-^°) (6) 

in the critical region. For brevity, only in this subsection we have set r = r(l; 0) = \ — K\l K\ c {\\ 0). 
The crossover from the 3D to the 2D critical behavior, as R — > 0, can be described in terms of 
a direct generalization 2 -! 3 -^^ of the usual phenomenological scaling hypothesis valid for isotropic 
systems. Specifically, it is assumed that, for sufficiently small r and R, the scaling form of the 
singular part /(l; r, h, R) of the free energy in a field h embodies also the anisotropy parameter R 
as follows 

/(l; r, h, R) « T 2-«(i;0) F(/lT -/3(i;0)-7(i;0) ( Rr -^ (7) 

where F is a universal function. The exponents a(l;0) and /3(1;0) refer to the specific heat and 
the magnetization of the 2D Ising model. The quantity </>, called crossover exponent, is universal 
and must coincide ^ 32133134 with 7(1; 0), the exponent of the susceptibility of the 2D Ising model. 
Taking two derivatives with respect to h in eq.© one obtains that the susceptibility in zero field 
is given by 

X (1;K 1 ,R)^A^t~^X^(B^Rt^). (8) 

where (x) , called universal susceptibility crossover-scaling function is uniquely defined by choos- 
ing the normalization X^°\0) = dX dx ^ = 1. and are non-universal scale factors. Here 
and in what follows, a superscript zero is attached to all quantities related to the 0th moment of 
the correlation-function. 

The scaling forms eqs.([7|) and ([8]) provide an interpolation between the critical behaviors in 
2D (i.e. for R = 0) and in 3D, for small non- vanishing R, and thus can describe both of them. 
In particular, by the normalization of X^°\x), eq.([8|) is consistent with the 2D critical behavior 
eq.© of the susceptibility in 2D. On the other hand, the consistency of eq.® with the 3D critical 
behavior xi^K^R) w A<® [r(l; R)]'^ 1 '^ is achieved by assuming that X(°\x) is sin gular as 
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K\ — > i^i c (l; i?) and, for x in a vicinity of x, has the structure 

[ ' (1 - x/i)7(lifl) W 

with 

i = J B(°) J R[r i? (l;0)]- 7(1;0) (10) 

and tr(1; 0) = 1 - K lc (l; R)/K lc (l; 0). Due to the universality of X^°\x), also the constants X (0) 
and x are universal. 

For small positive R, the solution of eq. ljlOjl yields the reduced shift of the critical temperature 
of the 3D Ising system with anisotropy R from the critical temperature of its 2D limit, which has 
the following asymptotic behavior 

K lc (l ] 0)/K lc (l;R)-l~R 1 /' t > (11) 

with (ft = 7(1; 0). Therefore this important result is a simple consequence of the crossover-scaling 
ansatz eq.flU) and of eq.(|5]), the critical behavior of x(l;ifi,0). 

The validity of the extended scaling assumptions eqs.jZJ and (jSJ) can be further tested by a 
numerical HT analysis of the asymptotic behavior as r — > of the successive partial derivatives 

4 j (i^i,o)-( — m — ) r=o (12) 

of x(lj-^i)-^) with respect to R, evaluated in the R = limit. The critical behavior of these 
quantities is defined by the asymptotic form 

H(°)(l;^ 1 ,0)«Ci°)(l)r-^ (13) 

as r — ► 0. By the extended scaling hypothesis eq.® the exponents of divergence X s should satisfy 
the relation 

A s = 7 (i;0) + ^ = (s + i) 7 (i;0). (14) 



For s = 1, eq. (|14p can also be seen as an immediate consequence of the relation 

E? ) (N;K l ,0) = 2K 1 [ X (N;K 1 ,0)} 2 (15) 



proved in Ref.pj] for N- vector models with arbitrary N. 

For N = 1 and s = 2,3, the validity of eq. (|14p is confirmed using the inequalities 

8/r?[x(l;tfi,0)] 3 > E< 0) (1;K U 0) > 4K 1 2 [ X (1; K u 0)] 3 (16) 

48^ 3 [ x (l;^i,0)] 4 > Sf (l;#i,0) > 8^ 3 [ x (l;^i,0)] 4 . (17) 

A proof 7 of eqs. (|16p and (|17p is based on classical correlation inequalities known to hold in the 
N = 1 case. Some generalization of eqs. (fT6|) and (fT7|) might still be valid also for models with 
N > 1. 

For the 1th moment of the correlation-function, one can assume the validity of the extended 
scaling form 

mW(l;Ki,R) « ^^(i^-Mi^m^)^^^)) (18) 

As a consequence, a generalization of eq. (|14p is thus obtained also for the exponents of divergence 
pL s of the successive i?-derivatives 

4 (1, ^ i, o) - ( ^ j H=o (19) 
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of the second moment of the correlation- function m^ 2 \l; Ki,R), which are defined by the asymp- 
totic behavior 

H( 2 )(l;#i,0) ^CP(1)t-^ (20) 

as r — > 0. From eq,(|18|) it follows that the exponents /i s should satisfy the relation 

Us = 2i/(l; 0) + (s + 1) 7 (1; 0). (21) 

For s = 1, the validity of eq. ([2Tj) is immediately proved^ by using an analogue of eq. (fT5|) for the 
second moment mS 2 ' (N; K i , R) of the correlation- function: 

E^iN-^^O) =2K 1 {[ X (N;K 1 ,0)} 2 + 2 X (N;K 1 ,0)m^(N;K 1 ,0)}. (22) 

Notice that, as for eq. (|15|) . also the validity of eq,(|22|) is not limited to the N = 1 model. 

For s = 2 and 3, inequalities analogous to eq,(|16|) and eq. (|17|) can be derived^ also for 
si 2) (l;ifi,0) 3 thus justifying eq.((3D]). A numerical test£±&±± of eas.([T3" ]) .(fI3 j) . (1201) and ([21]) gave 
support to the crossover-scaling ansatz for N = 1. 

As a final remark, let us point out that, using the variables Ki and R = 1/R, more convenient 
in the R — > oo limit in which the system becomes an array of one-dimensional spin chains, also the 
following relations, valid for arbitrary N, are obtained: 

■ amm(N;K 2 ,R) \ _ = { 2 + )m< s >(iV; Ki,0)} (24) 

ait /R=0 

Here x(iV; .K2, 0) and m^ 2 '(N; K2, 0) indicate, respectively, the susceptibility and the second mo- 
ment of the correlation-function of the anisotropic N- vector model for R = 0, i.e. in ID. Eqs. (|15p . 
(|22|) . (|23p and (|24|) are quite helpful also to validate the computation of the bivariate series expan- 
sion. 

It is now also clear— ^ how to compute an expansion of the universal susceptibility crossover- 
scaling function X^°\x) in powers of x. Observing that the critical amplitudes Cg°\l) in eq. lfTHj) 
are expressed in terms of (x) as 



(25) 

dx s / x=o 

and that the dependence on the non-universal quantities and disappears from the ratios 

r (o) r (0) 

* = >f • (26) 

the expansion of (x) for small x can be written in the form 

x(°)( x ) = 1 + x + —x 2 + 9l9i x 3 + 9lQMk x 4 + 9l9l9l9± x 5 + 9l9^9l9l9i x 6 ... (27) 

which is universal i.e. the coefficients are independent of the lattice structure. The calculation 
of (x) can be numerically extended to the whole interval [0, x] by Pade approximants. This 
approximation can be shown to yield an asymptotic description of the properties of x(lj-Kij-R) 
which, in the range of validity of the crossover-scaling ansatz, is consistent with those obtained 
from other approaches. 
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B. The N > 2 models 



A slightly different formulation of the extended phenomenological scaling is necessary in the 
N > 2 cases that will be considered in this subsection, simply because the asymptotic relations 
eqs. ([7]) , ([5]), and therefore (|13|) and (|2U|) cannot be valid anymore. Before introducing this 

issue, it is convenient to make a brief digression to recall how the critical behavior in 2D has 
been characterized^ in the N = 2 case. As K\ — > K\ c (2\§) from below, the divergence of the 
correlation-length of the 2D XY model is dominated by an exponential singularity 

£(2; #1,0) «U(2;^i,0) = Dexp[bT-°}. (28) 

In this and the following sections, for brevity, we have set r = r(2;0) = 1 — K\/ K\ c {2] 0). The 
universal exponent a is expected to take the value a = 1/2, while b is a non-universal positive 
constant. The critical behavior of the singular part of the free energy is predicted to be 

/ sins (2; #!,()) «/„(2;jri,0) = F[C as (2; K^O)}- 2 (29) 

with F a non-universal amplitude, while for the susceptibility one has 

X (2; K x , 0) « X as(2; 0) = A® r"^ exp [(2 - r/(2; 0))^]. (30) 

as i^i — > Ki c (2;0) from below. For > i^i c (2;0) both £ and x are infinite. The quantity 
7/(2; 0) = 1/4 is the exponent characterizing the large-distance behavior at criticality of the spin-spin 
correlation- function for the 2D XY model. In eq. (|30p . the presence of a multiplicative correction 
to the leading singular behavior by a power of the logarithm of £ (equivalently by a power of t~°") 
and the value of the exponent 9 are still controversial. A rediscussion— of the renormalization 
group approach indicates that 9 = 0, while a recent high-order HT study 30 and a high-precision 
MonteCarlo study 3 ^ support the estimate 9 ~ 1/16. However, the precise value of 9 is practically 
irrelevant in our discussion of scaling. The numerical value of the susceptibility critical amplitude 
A^ ) depends on the value assumed for 9, but it is also irrelevant in the determination of the 
crossover-scaling function X^°'(x). 

Let us now return to the scaling issue and make the natural assumption that the crossover- 
scaling ansatz, introduced in the Ising case for the singular part of the free energy in an external 
field of modulus h, can be simply generalized to the XY model case as follows 

f(2;r,h,R) « [£ as (2; K u 0)]- 2 F(hU{^ K u 0)[ X a S (2; K 1} 0)] 1 / 2 , R Xas (2; K u 0)) (31) 

for sufficiently small positive R and r. Taking two derivatives with respect to the field, we get the 
generalized scaling form for the susceptibility in zero field 

x(2;K 1 ,R)^Xas(2;K 1 ,0)X^^Rxas(2;K 1 ,0^ (32) 
Here (x) is a universal crossover-scaling function, that can be uniquely defined assuming that 

x®(o) = dX 2 {0) =1- # (0) 

is a non- universal scale factor. 
The generalized crossover-scaling forms eqs. (|3ip and (|32p are immediately shown to reduce to 
eqs. ([7]) and flSJ) in the N = 1 case, by using eqs.© and © and the scaling laws. However, for 
N = 2, these forms have the additional virtue of correctly allowing for the fact that the critical 
singularities in 2D are not power-like, since the exponents 7(2; R) and v(2; R) (as well as a(2; R) 
and /3(2; R)) are ill-defined in the R — > limit and also no crossover exponent (ft exists. Of course, 
the scaling forms eq. (|3ip and (|32[) can as well be written exclusively in terms of £ thanks to eq. 

HMD- 
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By the normalization of X^°\x), the scaling form eq. (|32p is consistent with the 2D critical 
behavior, since for R = and r — > 0, one has x(2;-Ki,0) pa Xas(2; K\, 0). On the other hand, 
for small but non- vanishing R, the conventional 3-D critical behavior %(2;ifi,i?) ~ ^4r~ 7 ( 2; ^ as 
r(2; i?) — > 0, is recovered by assuming that (x) has the same singularity structure as in eq.@, 
i.e. that X (0) ( x) ~ ^_ a . /j.\ 7 (2;.r) ; when x is in a neighborhood of x, with 

x = J B(°) J Rxa S (2;^i c (2;i?),0). (33) 

As in the N = 1 case, by solving eq. (|33j) . we can obtain the small- R asymptotic behavior of the 
reduced critical-temperature shift of the anisotropic 3-D model from its 2D limiting value 

K lc (2; 0)/K lc (2; R) - 1 » l//[ln(i?/VF)] 2 . (34) 

Here V = [2 - ??(2;0)] 2 6 2 and W = x/BA<® are non-universal constants. The higher-order 
corrections to the leading behavior in eq. (j34j) . are also expressed in terms of inverse powers of 
[lni?| and possibly of ln|lni?|, depending on the value of the exponent 9. It is interesting to point 
out that the original argumen t 3 l3S for eq. (|34p relied on approximate renormalization group ideas, 
rather than being a simple consequence of the generalization of the crossover-scaling ansatz and of 
the exponentially-singular critical behavior of the XY model. 

As a further immediate consequence of our scaling assumption eq. (|32p . the divergence of the 
successive ^-derivatives of the susceptibility turns out to have the structure 

E^(2;K U 0) ^Ci°\2)[ X as(2;K 1 ,0)] S+1 ^ r e ^exp[(2 - V (2;0))b(s + l)r^] (35) 

as r — > 0, which can be seen as a natural generalization of eq. (|13p . For s = 1, the validity of 
eq. (|35|) is an obvious consequence of eq. (|15|) . For s = 2 and 3, it might follow from some extension 
of the inequalities eq. (fl~6|) and eq. (fT7|) . Eq. ([35|) will be numerically tested in subsection IV C, by 
studying the behavior of si°^(2; Ki,0) as K\ — > K\ c (2\ 0), for the first six values of s. 

By assuming for the second moment (2; K\, R) of the correlation-function the crossover- 
scaling form 

m< 2 > (2; K u R) « m^{2;K u 0)X^ ^R Xas (2;K 1 ,0)) (36) 

where mi) {2; K ly 0) = 4[£ as (2; K 1 , 0)} 2 x as (2; K x , 0), the analogue of eq.jH]) is obtained 

3?) (2; K ± , 0) « (2) [U (2;K l , 0)} 2 [ Xa s {2;K U 0)] s+1 (37) 

By eq. ()22p . this equation is certainly valid for s = 1 and we shall suppose that it is true also for 
s > 1. Also eq. (|37p will be tested numerically in subsection IV C, by the same method used for 
eq.fSD. 

The HT expansions of ^^(2; K\, 0) and 3g (2; K±, 0) can be immediately read from our tables 
of the series coefficients of %(2; K±, K2) and m^ 2 ^(2; K±, K2), respectively. 

The expansion in powers of x of the universal scaling function X^ (x) can be computed by the 
same procedure already outlined for the N = 1 case. 

A generalized scaling assumption can be made also to describe the crossover of the 3-D system 
to a set of one-dimensional non-interacting XY chains as R — * 00. It is now convenient to shift 
to the variables K2 and R = 1/R. One has simply to observe that in this case K2 C (2; R) — * 00 as 
R — > and that, at criticality, the divergence of the susceptibility in one dimension is 

X (2;K 2 ,0) w Xas (2;K 2 ,0) = A^K 2 2 (38) 
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Then it is natural to introduce a universal susceptibility scaling function X^°\x), normalized like 
(x) and with the same singularity structure and assume that 

X (2; K 2 , R) « X as(2; K 2 ,0)X<® (S®Rxa.(2; K 2 , 0)) . (39) 

It follows that 

K 2c {2;R) -FT 1 ' 2 (40) 

as R — > 0. Moreover, predictions for the fi-derivatives of the susceptibility, similar to those 
mentioned above for the -R-derivatives, are easily obtained along the same lines. 

Let us finally point out that generalized crossover-scaling assumptions of the same form as 
eqs, ([31|) and (|32|) can be written down also for the AT- vector spin models with N > 2. The 
main difference is that for these models, K\ C {N\R) — > oo as R — > and that, in 2D, the critical 
divergence of the correlation- function moments and of the correlation- length is exponential in K\ . 
As a result, we can conclude that for small R, we have K\ C (N; R) ~ |lni?|, by the same scaling 
arguments used above. Here, we shall not further investigate the models with N > 2, but only 
note that a general discussion 39 of the numerical difficulties of the HT-expansion approach to the 
A^-vector models in 2D suggests that a HT study of the crossover might also meet with similar 
problems caused by the unphysical singularities in the complex inverse-temperature plane revealed 
by a large- N study^. 

IV. A NUMERICAL ANALYSIS 

We shall now turn to an analysis of our HT expansions to study the crossover behavior (|34p of the 
critical inverse-temperature K\ c {2; R) as R — > or R — > oo, to check the universality with respect 
to R of the critical exponents of the anisotropic system, to test the validity of the consequences 
eqs. (jMJ), ([3"5]) . ([3"7]) and ([30]) of the crossover-scaling assumptions eqs.([3"2~]). ([50]) . ([3"9"]) and finally 
to obtain an approximation of the susceptibility scaling function X^\x). 

A. Estimates of the critical temperature as function of R 

Let us first study the behavior of the susceptibility %(2; K\,R) as function of K\ at fixed values 
of R, to determine the critical locus. We shall later use these results to bias the computation of 
the critical exponents of the susceptibility and of the correlation-length and verify that they satisfy 
the universality hypothesis, as long as R > 0. For simplicity, we have analyzed the behavior of 
our expansions as functions of K\ at fixed values of R (or as functions of K 2 at fixed values of R) , 
using the conventional single-variable methods^ of series analysis, namely PAs or inhomogeneous 
differential approximants (DAs). It may be helpful to recall that in the DA approach a (single- 
variable) power series is resummed by expressing it as the solution of a linear (first- or higher-order) 
differential equation with polynomial coefficients and inhomogeneous term, appropriately defined 
in terms of the coefficients of the given series. We believe that taking advantage of the large 
number of series coefficients presently available, also more complex methods of series analysis, like 
multivariate PAs-^ or partial-differential approximants^, become now worth exploring. However, 
we have not yet thoroughly pursued these approaches. 

In our computation of the critical temperature, for R > 0.05, we have used second-order DAs 
and have considered the class of [k, I, m; n] DAs restricted by the conditions: 13 < k+l+m+n < 19 
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with k > 3; I > 3; m > 3. Among these, we have selected the DAs with the additional properties of 
being defect-free, i.e. having sufficiently isolated physical singularities, and of being near-diagonal, 
i.e. such that \k — l\, \l — m\ < 2, 1 < n < 4. Finally, we have not "biased" the DAs, i.e. we 
have not discarded the approximants yielding critical exponents with values outside appropriate 
limits. Here we are using the standard notation in which k,l,m,n denote the degrees of the 
polynomial coefficients and of the inhomogeneus term of the differential equation defining the DA. 
These rather technical specifications are given only to make our results completely reproducible. 
However, we have always made sure that our final estimates, within a fraction of their uncertainties, 
are essentially independent of the precise definition of the DA class examined. At a given value 
of R, our central estimate of Ki c {2; R) is the sample mean of the locations of the singularities in 
this class of DAs, taken after dropping evident outliers. A small multiple of the spread of the 
reduced sample is taken as an estimate of the uncertainty. Possible residual unsaturated trends of 
the central estimates, as the number of series coefficients used in the calculation is increased, have 
been accounted for by attaching generous error bars to the final estimates. Sometimes we have 
made small corrections of these central estimates (always well within the uncertainties) suggested 
by a comparison with the sequence of Zinn-Justin modified-rati o 16142 estimates. 

Coming to the numerical results, let us first note that, for R = 1, our estimate Ki c {2; 1) = 
0.22710(3) compares well both with our older estimate 28 K lc (2; 1) = 0.227095(10) obtained from 
HT expansions of order 21 for the isotropic 3D XY system and with the much more precise^ 3 - recent 
determination Ki c (2;l) = 0.2270827(5), obtained from a high-accuracy MonteCarlo simulation. 

The determination of the line of critical points Ki c (2; R) for very small values of R, as required 
for a numerical test of the predicted crossover behavior eq. lfMj) . is a delicate task both by series 
study and by simulation. The progressive decoupling of the horizontal layers as R — ► 0, makes 
a reliable extrapolation of the susceptibility to its genuinely 3D critical behavior possible only 
from a closer and closer vicinity of the critical point. (See also the comments to FigJS] in the 
next subsection.) As a consequence, the precision of the estimates of the location of the critical 
point and of the critical exponents tends to deteriorate as R — > 0. We have observed that, for 
R < 0.05, a reasonable increasing behavior of K\ c (2\ R) is still obtained simply by assuming what 
will be actually borne out by our analysis, namely that the critical exponent of the susceptibility 
is universal, to a good approximation, all along the critical locus for R > and thus by imposing 
also a "weak" bias on this exponent when computing K\ C (2]R). This simply amounts to discard 
from the sample of our data the critical temperature estimates obtained from DAs whose exponent 
at the singularity differs by more than 10% from the expected value of 7(2; 1). Even after this 
simple improvement of the analysis, for R < R mm ~ 0.0015 , our expansions do not anymore 
seem to be long enough to locate the critical point with acceptable accuracy. Therefore we shall 
not report estimates for R < R m %n- Our results for the reduced shift S = Ki c (2;0)/Ki c (2;R) — 1 
of the critical temperature from its 2D limiting value are plotted vs R 2 ^ 3 in FigfTJ The value 
Ki c (2;0) = 0.56000(5) of the critical inverse-temperature of the 2D XY model on the square 
lattice has been taken from recen t 29 ' 30 ' 36 high-precision studies. 

A short list of our numerical DA estimates of Ki c (2;R) for 0.005 < R < 3.4 can be found in 
Table E 

It should be noticed that the critical curve must be separately symmetric under the transforma- 
tions K\ — ► —K\ and K2 — ► —Ki. Therefore the ferromagnetic phase diagram, so far represented, 
should be completed by the remaining branches of the critical locus. 

In Fig(T]we have also shown that S has a very simple dependence on R, valid to a high accuracy 
over a wide range of intermediate-large (but not too large) values of R. The expression f(R) = 
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aR 9 + c with a ~ 1.245, g ~ 0.661 and c ~ 0.221 interpolates quite accurately our data points 
for S in the interval 0.025 < R < 2.9, visibly departing from them only for small R, i.e. in the 
region where the crossover is expected to occur. The value of the exponent g justifies our choice to 
plot S vs i? 2 / 3 . It would be interesting to look for an analytic explanation of this purely empirical 
remark. 

It is interesting to remember that a not very different i? 4 / 7 law was conjectured 4 ^ to describe 
the i?-dependence of the critical temperature over the region 0.01 < R < 0.7, while for R > 0.7, a 
linear law was proposed. This suggestion was based on a self-consistent mean-field approximation, 
which is probably not very accurate for small R and is certainly inaccurate for intermediate-large 
R. 

Let us now turn to the small-R region to study the features of the crossover behavior. Figj2] 
is a blow-up of the lower left corner of FigfTJ showing our estimates of the critical temperature 
for R < 0.13. We can notice that they clearly depart from the simple i? 2 / 3 behavior indicated 
by a long-dashed line. A continuous line indicates the result of a fit of the asymptotic expression 
Vy[ln(it!/iy)] 2 of eq. ([3"4"j) predicted by the crossover-scaling theory to describe the behavior of the 
reduced temperature shift in the crossover region. Restricting the fit to the critical-temperature 
estimates which fall within the window (0.0025, R ma x = 0.1), we determine the values of the pa- 
rameters V ~ 11.34 and W ~ 12.7. In the same figure, for comparison, we have reported also some 
old MonteCarlo estimates 21 of K\ c (2\ R) in the small R range. In spite of their order-of-magnitude 
agreement with our series estimates, these simulation data seem to suggest a qualitatively different 
behavior as R — > 0. We must therefore suppose that they are affected by large errors (unfortu- 
nately not assessed in Ref . j2l| ) from underestimated finite-size effects, because the volume of the 
simulated system was small (at most 20 3 ) and no finite-size-scaling analysis was performed. A 
more recent simulation^ 2 ., also carried out with 20 3 sized systems, and extending to much smaller 
values of R, is likely to suffer from similar problems, although its results seem to show a better 
agreement with ours. 

Some remarks have to be made on Figj2j Firstly, one should be aware that a logarithmic 
behavior is quite difficult to identify numerically by using data which refer to only a two-decade 
variation of the independent variable, since, over a restricted range, it can also be well represented 
as a power-law behavior with a small exponent. Moreover the choice of the window of values of 
R to be studied, is a delicate issue, because the upper end of the window should be sufficiently 
small that both the crossover-scaling assumptions and the asymptotic behavior eq. (|34p apply with 
small corrections, while the lower end should be sufficiently large that our estimates of the critical 
temperature are not too uncertain. In our case the value of R ma x can be varied by a factor of two 
or more, still obtaining good fits of the same functional form, with not very different values of the 
parameters V and W. Finally, we observe that the best-fit value of V is not very different from 
its expected value (2 — ry(2; 0)) 2 6 2 w 9.5, while the value of W is much larger. It is reasonable to 
interpret these results as an indication that R m in is still too large and therefore the higher-order 
corrections to the asymptotic form eq. (|34p . suppressed only by inverse powers of |m-R|, are still 
important, so that V and W can only be effective parameters. To conclude, in spite of the notable 
extension of the HT expansions we have analyzed, these results still can give only a suggestive 
indication that our estimates of the critical temperature for small R are compatible with the 
predicted asymptotic behavior eq. ([34"|) . 

On the other hand, one might wish to describe also this crossover behavior by a power law and 
fit an expression f(R) = aR 9 ' , of the same form as eq. (|1 lj) valid in the N = 1 case, to the same 
small- R sample of our data points. We then find the following values of the parameters: a ~ 1.08 
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and g' ~ 0.354. We have displayed also this fit in Figj2j Thus, if for consistency we describe also 
the critical behavior of the 2D XY model in terms of conventional power laws, this result suggests 
an unusually large susceptibility exponent, i.e. 7(2; 0) = 1/g' « 2.8. We can also observe that, 
subdividing the small-i? range into two intervals, this kind of fit would yield a smaller exponent 
g' in the leftmost interval. This suggests that even smaller values of the exponent g' (and thus 
larger values of the susceptibility exponent) are likely to be found, if it were possible to further 
reduce both ends of the window of values of R under consideration, for example by using future 
significantly longer HT expansions. It is then reasonable to believe that this power-law crossover 
is only apparent, because the exponent is independent and vanishes as R — > 0, and hence that we 
are actually representing a genuine logarithmic behavior as a power-law behavior. 

In Figl3l we have plotted {2K 2c (2; R))^ 1 vs R 2 ^ 3 . The behavior of the curve as R — > also gives 
some support to the validity of eq. (|40p and therefore it confirms the generalized scaling approach 
to the crossover from 3D to ID. 



B. Universality of the critical exponents with respect to R 

Unfortunately, so far we have computed expansions only for the sc lattice structure and therefore 
our verification of the universality of certain quantities does not include their lattice independence, 
but is limited to a test of their independence on R. 

In the range R > 0.05, in which our computation of the critical temperature was completely un- 
biased, we have also evaluated both the critical exponents j(2; R) and v(2\ R) by using second-order 
DAs, chosen in the above specified class and biased with our estimates of K\ c {2\ R). We have also 
determined, by simple first-order DAs, the ratio of the logarithmic derivatives of m^ 2 \N; K\, R)/K\ 
and x(iV; K\, R), whose value at K\ C {2;R) yields the ratio u(2; R)/j(2; R). In FigHJ we have 
plotted vs R our estimates of these exponents and of their ratio normalized to our estimates of 
7(2; 1) = 1.328(10), i/(2; 1) = 0.679(8) and i/(2; 1)/t(2; 1) = 0.5099(6), respectively. The central 
values of our estimates of the normalized exponents and of the exponent ratio are very near to unity 
and fairly independent of R within ~ 0.5%, over a range wider than that shown in the figure, except 
for very small or large R, due to the expected crossover. These results support the universality with 
respect to R of the critical behavior for the anisotropic 3D model. It is fair to note that, while the 
deviations of our exponent estimates for the anisotropic system from those for the R = 1 system 
are uniformly very small, our estimates of 7(2; 1) and z/(2; 1) are larger, by < 1%, than the most 
recent 4 ^ estimates 7(2; 1) = 1.3178(2) and u(2; 1) = 0.67155(27). On the contrary, the deviation 
of our estimates of the exponent ratio from its recent determination z/(2; l)/7(2; 1) = 0.5096(3) 
is quite small. Substantially longer series or improved methods of exponent determinatio n 26 ' 42 ' 45 
might be necessary to remove the discrepancy in the exponent estimates, which certainly can be 
ascribed to the residual influence of the corrections to scaling. 

In Figj5]we have plotted, for various fixed values of R, the effective exponen&^L ^ e ff(2;K\,R) 
of the susceptibility 

dln X (2-K u R) 

leff (2, Kl ,R)-- dinr{2 _ R) . (41) 

Roughly speaking, this quantity represents the local value of the critical exponent which would 
be inferred by a measure of the susceptibility in a neighborhood of K\ and only in the critical 
limit it coincides with the asymptotic critical exponent. By showing how the local critical behavior 
changes, this quantity is helpful to describe the crossover phenomena. The curves shown in the 
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figure are obtained simply by forming the highest-order, defect-free, diagonal or near-diagonal PAs 
of the HT expansion of the right-hand side of eq. (|4ip . The figure illustrates the narrowing down of 
the genuinely 3D critical region as R — > 0, that was already discussed in the preceding subsection. 
More precisely, when R is not too small, the effective exponent is approximately independent of the 
temperature and stays close to its 3D asymptotic value. On the other hand, going to sufficiently 
small R, the effective exponent becomes increasingly temperature-sensitive: in most of the HT 
temperature range it takes rather large values, which somehow reflect the exponential 2D critical 
singularity, and only as K\ gets very close to the critical value, it approaches the 3D asymptotic 
value. 

In Fig JS] we have plotted, for various fixed values of R, the effective exponent r y e ff(2;K2,R) 
of the susceptibility, defined in strict analogy with eq. (|4ip to describe the crossover from 3-D to 
ID. Also in this case, the curves indicate a transition from a region of very high values of the 
effective exponent, reflecting the critical divergence of the susceptibility of the ID system, to the 
asymptotic region where the 3D critical behavior is attained. 



C. Critical behavior of Eq (2; K\, 0)). An approximation of the susceptibility 
crossover-scaling function X(°\x). 

As sensitive and specific indicators to test the critical behavior eq. (|35|) of the successive R- 
derivatives of the susceptibility (2; -K"i, 0), we have formed the normalized ratios of their loga- 
rithmic derivatives 

Gf _ _J_ rfl n[^»(2; Kl ,0)]/ (igl 
o + ldl„[4°»(2;A' 1 ,0)]/ t iA' 1 

If eq. ([35|) . showing the critical behavior of (2; K\, 0), is valid, the quantities G s (Ki) are 
expected to tend to unity, as r — > 0, independently of s. This result is obvious for s = 1, thanks 
to the exact relation eq. (fl"5j) . but for s > 1 it ought to be numerically tested. We have used 
near-diagonal first-order DAs of the HT expansions of G^(Ki) to estimate the limits of these 
quantities as K\ — > _Ki c (2;0). Our results, summarized in the Table |H] support the validity of 
the asymptotic relation eq. (|35p to a good approximation when s = 2,3,4,5,6 and therefore lend 
support to the validity of the generalized crossover-scaling assumption eq. (|32p . In general, it should 
not be surprising that the precision of these, as well as of the following estimates, is smaller than 
that reported in previous studies of much shorter series in the N = 1 case, in which the critical 
singularities are power-like, because of the more complex structure of the critical singularity in the 
2D XY model. One should also notice that the number of nontrivial expansion coefficients of the 
.R-derivatives of the moments of the correlation-function and hence the precision of the numerical 
approximation for the quantities related to them, decreases as s increases. Therefore these tests 
are less reliable for values of s larger than those examined here. 

A completely parallel study of the indicator Gs (Ki), a strict analogue of Gs{K\) for the 

quantity 7g/ 2 -Xi oj] s ' which can be associated to the i?-derivatives of the second moment of the 
correlation- function, also leads to results in agreement with eq. ([37|) . albeit slightly less precise, 
because of the known slower convergence of the second moment expansion. Thus also the validity 
of eq. (f36"l) is confirmed. Our estimates of the limits of the quantities G s (Ki) as K\ — > iCi c (2;0) 
are also reported in the Table ITI1 

Let us now turn to the study of the universal crossover-scaling function (x) of the susceptibil- 
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ity, following the procedure^ described above in detail for the N = 1 case at the end of subsection 
III A. As a first step, we have to determine the critical amplitudes ci°^(2) of 3g (2; Ki,0), defined 
in eq.422) by the limit of the effective amplitudes d° ) (2;K 1 ,0) = E { ° ) (2;K 1 ,0)/x s (2;K 1 ,0) as 
K\ — > K\ c {2; 0). The amplitudes C^°\2) are then used to form the universal ratios Q s , defined by 
eq , ([26]) . Alternatively, the ratios Q s can also be determined, directly and perhaps more accurately, 
by extrapolating the HT expansions of the effective ratios 

Q. m ^« K ^«™ (43) 

to K\ = Ki c (2;0). The estimates so obtained for the universal quantities Q s = 
'&aijc 1 -*K-L C Qs(^]Ki,0) are reported in Table UTTl 

The coefficients of the small- x expansion of A(°) (x) are finally expressed in terms of the Q s , as 
indicated in eq. lffi 



X i0 \x) = l + x + 0.792(3)x 2 +0.600(5)x 3 +0.44(l)x 4 + 0.32(l)x 5 +0.23(2)rr 6 + 0.16(3)x 7 + ... (44) 

Having assumed that X^°\x) is singular at x and has the form X(°\x) ks WZ^T^W^) m a vicinity 
of x, we can give a reasonably good estimate of x by locating the nearest pole of the highest-order, 
defect-free, near-diagonal PAs of the expansion of [X(°)(x)] 1 /7(2;-R). Choosing j(2;R) = 1.3178(2) 
and allowing for the uncertainties of the expansion coefficients in eq. (|44p and the spread of the 
PA singularities, we can estimate x = 1.475(15) and X^ = 1.154(15). A direct estimate of x, by 
extrapolating to R = the argument BRx a s(2; K\ c (2\ R), 0) of the scaling function, fails to yield 
a more accurate result because of the extrapolation uncertainties. 

Let us represent the scaling function simply as X^°'(x) = P{x/x){l — x/x)~'~ 1 ^ 2 ' ,Ft '\ where P(x/x) 
is some function interpolating between the small x and the large x behavior of X^ (x) and therefore 
taking the values P (0) = 1 and P(l) = X^°\ It turns out that the simple linear expression 
P{x/x) = 1 + (A(°) - l)x/x is a quite accurate approximation of the regularized scaling function 
X(°\x)(l — x / x)^( 2 ' R \ We have then used this approximate form of the crossover-scaling function 
to compute the effective exponent 

K lc (2;R) t(2;R) P'(z) j{2;R) ^ 

7e ^ (2 ' R) = g lc (2;0) [r(2;0)]^ ba[1 + Z{ W) + T^ )] (45) 

with z = x/x. The effective exponent computed using eq. ([4"5j) is plotted vs r(2; R) in FigfTJ in which 
we have shown the curves corresponding to the eight smallest values of R chosen in Figj5j The 
agreement between Figj5] and FigfT] can be considered as another satisfactory confirmation of the 
validity of the crossover-scaling ansatz, if one observes that the former is a good PA representation 
of the effective exponent on the whole interval < K\ < K\ c (2\ R), whereas the curves of FigfTJcan 
be quantitatively reliable only in the small range of validity of the crossover-scaling form eq.(|32p. 
for example when r(2; R) < 0.05 and R < 0.05. 



V. CONCLUSIONS 

We have added 187 coefficients to the already known 66 coefficients of the bivariate HT ex- 
pansion for the spin-spin correlation-function of the 3D XY model, with directionally anisotropic 
couplings, on the sc lattice. Analyzing these data by PA and DA methods, we have determined to 
a good precision the critical locus of the system in the ferromagnetic region and checked to a fair 
accuracy the universality of the critical exponents of the susceptibility and the correlation-length 
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with respect to the anisotropy parameter R. We have also shown that the main predictions of the 
extended scaling theory for the crossover from the 3-D to the 2D critical regime, concerning both 
the behavior of the line of critical points K\ c (2\ R) in the limit of small R and the critical divergence 
of the successive A-derivatives, at R = 0, of the susceptibility and of the second moment of the 
correlation- function, are compatible with the numerical extrapolations of our extended expansions. 
Finally, combining the HT expansions with the crossover-scaling ansatz, we have obtained a con- 
crete approximate representation of the crossover-scaling function X^°\x) for the susceptibility, 
and have shown that it reproduces, in an appropriately restricted temperature range, the effective 
exponent as computed by PAs on the HT side of the critical point. 



VI. APPENDIX 



For the Hamiltonian H an of eq. (pQ) ( anisotropic nn model) with z-anisotropic nn interactions in 

3-D, the HT expansion of the energy (and of the specific heat) can be simply formed in terms of 

the nn correlation in the direction of the z-axis: 

1 1 479 321 11 

C(000, 001; 2; K x , K 2 ) = K 2 +4K 2 1 K 2 --Kl+32KtK 2 +-Kl+—KfK 2 +—KfKl-K 2 Kl--K 7 2 + 

^KfKl - ™°KiKi - ±KM> + —•<? + 503352137 K?K, + ^996529 

27 12 36 1 2 45 1 2 3024 2 6720 1 5760 1 2 

3577027019 8389362923 g? 7+ 3835391 g . 8+ 83^ g? u j57_ 2 13 _i0136?_ 

576 1 2 6912 1 2 1728 1 2 1920 1 2 8640 1 2 1935360 2 

77461196851 335838080671 , 4 » 6966586038139 , 2 . 113783648353 , n - 

241920 1 15120 1 2 77760 1 2 2592 1 2 

20055410369 _ 6476551 _ 442849 _ ^483_ 946523 17 

10368 1 2 1440 1 2 12960 1 2 90720 1 2 26127360 2 

87581038655119 M 2276495337319411 16 3 656734833787141 14 5 275625926210833 12 7 

Ai Ko-\ Ai A 9 H A-. A 9 H Ai An + 

65318400 1 14515200 1 2 604800 1 2 259200 1 2 

196795509446797 109 63925399531 811 _ 242287049 613 144905431 415 

1296000 1 2 57600 1 2 259200 1 2 3628800 1 2 

45679 r ,, r , 17 65467219 - q 3182306618881 ~ n r ^ 34008570927861617 , s - 
-A? A 9 17 A, 19 H A 20 A 2 H A| 8 A|+ 



7257600 1 z 2612736000 z 576000 1 32659200 

84027770340130537 r , 1fir ,= 13557570274734337 , 4 - 218753264133841 r q 

A^A^ H A| 4 A 7 H K\ 2 KZ+ 

7257600 1 2 680400 1 2 34560 1 2 

39683430286691 1011 _ 114458157007 813 4455945869 615 _ 165611879 417 

144000 1 2 259200 1 2 907200 1 2 2419200 1 2 

41921 KlKf + 249045899 ^ 21 + 



32659200 1 z 14370048000 

and of the nn correlation in the direction of the x-axis: 

3 1 31 479 321 1 

C(000, 100; 2; K u K 2 ) = A 1+ -A 1 3 +2A 1 A 2 + -A 1 5 +32A 1 3 A 2 - -K 7 l+ —KIk 2 + —K\k\- -K x K%+ 

2 6 48 2 4 b 

Z^J^+1631tfM+2162tfM+4^ 
120 1 12 1 2 1 2 12 1 2 1440 1 72 1 2 36 1 2 
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15389 9 
27 



KM _ _ * KlK n _ ^1^3 + ^ A -n A .| + 31037915 

*i 1 2 72 1 2 360 2 5040 1 10 1 2 72 1 2 

6002606 76 50104 5 8 1169 3 10 1 12 75180487 15 503352137 132 

27 ^1-^2 g ^1-^2 180~ 1 2 " 540 ^ 2 " 645120 Al + 1920 Al ^ 2 + 
16776996529 17885135095 8389362923 ^35391 5 83969^i 2+ 

3840 1 2 3456 1 2 13824 1 2 5760 1 2 11520 1 2 

457 r ^ 14 6506950039 , 7 77461196851 15 2 335838080671 13 tM 6966586038139 _. n 
120960 2 26127360 1 60480 1 2 8640 1 2 77760 1 

568918241765 9 8 20055410369 7 ^ 10 6476551 442849 3 14 2483 16 

20736 Al ^ 2 + 25920 Al ^ 2 5760 ^ 2 ~ W720 ^ ~ 1451520 Kl * 2 + 
-1102473407093 iq 87581038655119 17 2 2276495337319411 ^ 4 656734833787141 ,„ 
2612736000 1 14515200 1 2 7257600 1 2 518400 1 2 

275625926210833 ,, » 196795509446797 ~ rXin 63925399531 -, 12 242287049 . , 4 

K'E H iffif, H K7i^ 2 + 

345600 1 2 2592000 1 2 172800 1 2 1209600 1 2 

144905431 3 16 45679 18 _ 6986191770643 21 3182306618881 192 

29030400 1 2 130636800 1 2 14370048000 1 115200 1 2 

34008570927861617 17 4 84027770340130537 fi 13557570274734337 ,„ s 

14515200 1 2 5443200 1 2 777600 1 2 

218753264133841 n 10 39683430286691 g 12 114458157007 7 14 4455945869 5 16 

57600 1 2 345600 1 2 907200 1 2 4838400 1 2 

-165011879 41921 

21772800 1 2 653184000 2 

The HT expansion of the susceptibility x(2; K\,K-i) reads: 

X (2; Ki, Jf 2 ) = 1 + + 2if 2 + 12if 2 + 16JfiJf 2 + 2X 2 + 34Jf? + 80KfK 2 + 32X 1 X 2 + Jf|+ 

+ 328KfK 2 + 2^KlKl + AQKiKl + ^-Kf + \\MK\K 2 + XZTlKlKl + 468if 2 iff + ?>2K X K%+ 

--K$ + 529if? + HI5?jf6jf 2 + 6416^^1 + 3656X?X 2 3 + ^K\k\ + —K Y K\ - -K%+ 
3 3 3 6 

14933 ^ 7 + 36484^,6^ + 26838K 5 R 2 + 2 2\%\K\kI + 6985if 3 if 2 4 + %TlK\Kl - ^K X K Q 2 + 

1 T ^ 7 5737 o 107813 ^ fi ^ 2 344200 . ■> „_ „^ 4l .^ 4 30092 - 

— if J H Kf H ^7^2 + 102070^?^ H i^if 3 + 54480.K?!rf H K^K%+ 

24 2 2 1 3 1 12 3 12 1 2 3 1 2 

1202 23 „ „ 7 1 o 389393 ~ 727969 1585027 , 

if 2 iff ^i^2 + — H K? + 102094i^f^ 2 H H iff if? + 

3 12 3 2 12 2 60 1 1 2 1 2 3 1 2 

349093if 1 5 if| + ^I^iff + ™if?if 2 6 + ^fif 2 if 2 7 - 4X^ 2 8 + 

2608499^ 10 + 421|53^^ + 3675469^ + 6689645^^3 + 3850939^+ 

3 1 2 6 1 2 3 1 2 3 1 2 15 2 360 2 
3834323^n + ^I^o^ + 3555^32317^^ + 78955291^^ + 86092934^ + 

145884614 41039243 1607615 38584^ _ + ^^^o + 

27 1 2 27 1 2 9 1 2 9 1 2 90 1 2 45 2 

41 n 1254799 12 176710657 n 551185898 102 489144733 93 
"2160 X2 + ^^^ + 90 Kl K2+ 45 Xl ^ 2 + l^ Xl/f2 + 
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783263929 8 4 849262729 7 5 327546343 6 6 62191328 5 7 2922371 4 8 
^1^2 H ^ ^1^2 + 27 K i K 2 H ^ ^1^2 H — ^1^2 + 

364 793 2 10 283 n 1 12 84375807 13 226326302 

~"15~ 1 2 ~ "15~ 1 2 + 270^ 2 - 540 X2 + —560~ ^ + 45 Xl 
147264673^,^ + 6933895627^ + 3332454773^^ + 330361671^ | + 

17920524733 4843721705 H2?5939 1891861^ _ 3?7099^io + 

216 216 1 2 6 1 2 18 1 2 180 1 2 

499 2 n 58 191 3 6511729891 5315902873 

"20-^ 2 ~135^ 2 + 30240^ 2 + 20160 Xl + 420 ^ ^ 2+ 
4839028456 12 2 17700580409 n 3 44549564471 10 4 215942560781 9 5 

45 Kl * 2 + 45 Kl K * + 60 Xl K * + 270 + 

108374129725 86 19538343311 77 7458043925 68 766009373 B9 
216 12+ 1Q8 12+ 216 1 2 270 1 2 

5860253 _ 261878 5842 _ 2707 _J57_*u + 

180 1 2 135 1 2 135 1 2 3780 2 120960 2 
66498259799 r ^ 15 17559467153 , 4 ^ 1030585285427 ,o 2 1241088903283 , 2 - 
96768 1 560 1 3360 12 960 12 

545333177917^ n ^ 4 15628446835681 , n r ^ 11934250539301 9 . 4368379416899 o ~ 
192 1 2 4320 1 2 4320 1 2 3456 1 2 

127475570107 78 191978104529 69 3117156659 510 _ 188490343 411 
384 1 2 4320 1 2 1440 1 2 8640 1 2 

4246757 Hi™ _ ^57_ _ _J09_ 1B+ 

8640 1 2 6048 1 2 2520 2 193536 2 
1054178743699 16 1853158502071 15 5175865133419 14 2 31158713948009 13 3 

725760 Kl + 24192 Kl ^ + 6048 ^ ^ + 7560 Kl * 2 + 

2498158420849 , 2 , 6679846871285 „ n , 7582320713789 , n - 17084567411611 ~ - 
240 1 2 432 1 2 540 1 2 2160 1 2 

1162221448157 00 2235921456439 - ~ 6389841029 48134516 

432 K * Kl + 4320 KlKl + ^35— K ^ + ~ 45" + 

83823949 9312019 _ 233021 3257 „ _ _J483_ 

2160 1 2 15120 1 2 30240 1 2 17280 2 1451520 2 
63505993331 17 4789948744531 16 566876641827121 15 2 773547561026857 14 3 
13063680 Kl + 25920 Kl * 2 + 241920 Kl * 2 + 60480 Kl ^ 2 + 
21204512261149 ^ 1W4 4893241012078331 T ^ 12 T ^ 2601671974782203 „ n , 
576 1 2 77760 1 2 38880 1 2 

585636405603727 107 994750117740587 98 254490705579917 89 
12960 1 2 51840 1 2 51840 1 2 

17805962398507 „ 7 „m 261786374297 . 568937123 - tM 877704563 . , 3 

25920 1 2 6480 1 2 25920 1 2 36288 1 2 

27164849 3 14 _ 11452451 2 15 18229^ 16 _ 16151 17 19830277603399 18 
36288 1 2 725760 1 2 90720 1 2 26127360 2 3110400 1 

1439510708686567 , 7 2287661692966927 , 6 . 4687330747502107 , 5 o 

iVi Kn -\ ivi iv 9 H K-\ K -\- 

3265920 1 362880 1 2 120960 1 2 

3810140897792597 , 4 , 66815169002499599 T ^ 13 . 46842552495817169 T ^ 12 , 

K] K~ H K} 3 K% H K} 2 K5+ 

30240 1 2 272160 1 2 155520 1 2 

4658208957065401 , x , 6388084191664099 , o 6411755638079 ^ 9 „ 

Ki iv 9 H K-, K H iv!iv 9 + 

19440 1 2 51840 1 2 160 1 2 

133534055859457 810 1997340559319 711 261618690883 612 10212267019 513 

17280 KlK2 + 2592^ K ^ 2 + 10368^ K ^ 2 18144 K ^ + 



39863 
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49463021 4 14 30377351 3^15 111247 ^ 16 , 95791 ^ 17 45679 yl8 , 
90720 1 2 + 120960 1 2 17280 1 2 + 3265920 1 2 + 130636800 2 + 
8656980509809027 19 16990346289372257 18 181822163330740067 172 
653184000 1 + 16329600 1 2 + 10886400 1 2 + 

277976926650925547 16 3 253619834611062257 15 4 139504244168782903 14 5 

2419200 Kl * 2 + 604800 1 2 + 151200 Al ^ 2 + 

586373897113473881 6 154458497864612623 r ^ 12 7 188558667394136629 n r ^ 8 
453600 1 2 + 129600 1 2 + 259200 1 2 + 

188013182402455507 10 9 982091938311529 9 10 908320442109053 8 n 
648000 Kl K ' 2 + 13500 KlK2 + 86400 ^ 2 + 

30871722528221 712 6909844908217 613 131184772823 514 3786229351 415 
43200 1 2 + 907200 1 2 226800 1 2 + 302400 1 2 + 

167790391 3 16 18884699 2 17 _ 81803 18 529301 19 
806400 1 2 3628800 1 2 1166400 1 2 1306368000 2 
2985467351081077 - 397290645875080747 ^ 31654153295670463 , 8 - 

K-, H Ai Kl H Ai Kn-\- 

108864000 1 163296000 1 725760 1 2 

1091318923155197137 17 3 234970051697118173 16 4 18262470593869394021 15 5 

3265920 Kl K2 + 172800 Kl * 2 + 5443200 Kl K ' 2 + 

14461460147180084629 , 4 ~ 5092048288872136451 , , - 3104846020933512733 , 2 „ 

K} 4 KS H K\ z Kl -\ K} 2 K% + 

2721600 1 2 907200 1 2 777600 1 2 

1854733050821010329 n 9 96563772320850199 10 10 37397800587256709 9 n 

972000 1 2 162000 1 2 324000 1 2 

3185944243721879 o , 2 2088234805193 - , 3 1262744580547 , , 4 

259200 1 2 4032 1 2 226800 1 2 

222365630579 .,. 40181717827 . , fi 4507243337 , 7 

K?Kl 5 + KfK? 6 K?Ki 7 + 

907200 1 2 3628800 1 2 16329600 1 2 

174167243 _ 8229559 41921 

32659200 1 2 163296000 2 653184000 2 

811927408684296587^ 21 4385492759351981^ 20 ^ + 12241208160951859717^ 19 ^ 2 

14370048000 1 + 777600 1 2 108864000 1 2+ 

62360406854392578737 , 8 ~ 3128332273134190457 , 7 < 17192353387793409331 , fi . 

65318400 1 2 725760 1 2 1451520 1 2 

4073641693391130869 ~ 548535341435511268981 r „ 14r „ 7 149677215569420577263 „ 8 

193536 1 2 21772800 1 2 7257600 1 2 

22488614102184820471 r ^ 12 r ^ 9 7882602902572789 11 10 698300238344288257 n 

1944000 1 2 + 1800 1 2 + 648000 1 2 + 

721052835977723 „ 12 88211417051098193 o 253014314268331 Wr „ 14 

iff if, H if? iff H KJKo + 

4500 1 2 7257600 1 2 1036800 1 2 

-5964052451177 6 15 85859357495 5 16 11197132927 4 17 
580608 1 2 + 870912 ^ 2 + 4665600 1 2 + 

-5286752777 3 18 _ 39273503 219 36697 20 _ 168103 21 

65318400 1 2 46656000 1 2 13608000 1 2 1368576000 2 "' 
The expansion of the second moment of the correlation-function m( 2 )(2; Ki,K 2 ) reads: 

m< 2 )(2; K U K 2 ) = 4K ± + 2K 2 + 32if 2 + 32K 1 K 2 + 8K 2 + 162if 3 + 272if 2 if 2 + \28K X K 2 + 17iff+ 

yoyo 

Q72K 4 + 1680if 3 if 2 + \2A&K\kI + 336^i^| + 24K% + ~^-K\ + 85UKfK 2 + 

71 24772 114128 

9100if 3 if 2 2 + 3892X 2 if 3 + 640^1^ + — K\ + -^—Kf + — - — K\K 2 + 54848if 1 4 X 2 2 + 
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33552*3*5 + 9088*?*! + + f *« + 2™»*J + ^*»* 2+ 

287078^^1 + 236593^-fff + 93241^^ + 16622^^1 + ^^K^f + ^if 2 7 + 77996^+ 
HH^?^ a + ™^*f + 143 71041^f + + 

!^ A -= Al + ^*,** _ II*? + ™* ? + ^*«* 2+ 

U6|M5 A . 7A1 + 23325767^ + + 5953681^+ 

369267RtAl + + ™> KlK , _ ™«, + ™ K ;. + 

34012382 + 70813492 A ,, A .; + 3e3i3sg6K 7 K 3 + 33 223346*«X 2 4 + 

^*?*| + 4215454*;*. + ISf"*?*! + ^*?*| + ™ Kl A} + 

^187 611969977 985452181 1621421069 

90 2 360 1 45 1 18 1 2 

1577314567 T o, 1668258680 - . 3039023000^^, 1061738147 „ e 

9 1 2 9 1 2 27 1 2 27 1 2 

67469699 6031258 1^31267 _ 

9 1 2 9 1 2 90 1 2 45 2 

395 n 202640986 12 612509213 n 14744551642 102 
432 2 45 1 9 12+ 45 12 + 

11275330394 „ , 8562260260 8 . 19108224278 7 . 2830398088 „ B 
15 K\Kl + K\Ki + K\K\ + K\k\+ 

2196959680^ + 101314112^ + 3327988^ + ^^io + 

-^M 2 " + ^Kf + 58900571047 9232198298 
135 2 270 2 5040 1 45 1 

205863030521 112 549781840789 103 81925739897 94 73086198751 85 

180 12+ lg0 12 + lg 12+ lg 12 + 

480629562289 7 fi 160430759525 fi 7 2587900213 - 8 256823701 , „ 

216 + 216 ^ 2 + ~ LS^ ^* 2 + -^cS— 

91921193 _ 39139 _ 6856 22199 

180 1 2 12 1 2 135 2 30240 2 

3336209179 , 4 42155577203 , 3 r ^ 19258430024 12 2 177682250558 ^ n . 

112 1 70 5 1 2 15 1 2 

307643703989 2911760805313 „ - 768911042053 „ , 319585153487 ~ ~ 
15 + 7^ + U K\Kl + KlKl + 

27117476941 9822059771 135760609 6957800 3 n 

18 1 2 45 1 2 9 1 2 27 1 2 

837212 2 12 106661 13 _ _9959_^ 14 1721567587879 15 1740519488525 14 

135 1 2 1890 1 2 30240 2 23040 1 1008 1 2 

126704980674473 r ^ 1W2 42345007132019 , 2 - 84243175256177 < 

10080 1 2 960 1 2 960 1 2 

463449944770897 , . 361807221213917 n . 145911608744899 = T ^ 46852875522355 . . 

iVi iv 9 H a, iv 9 H iVi iv 9 H iVi iv 9 

4320 1 2 4320 1 2 3456 1 2 3456 1 2 

11447940964657 409797896227 . ,„ 111235514009 r , r ,, 35705321 ~ r , 2 

4320 + 1440 ^ 2 + 8640 ^ 2 + HiSSO"*^ ' + 
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122203853 213 134969 14 185291 15 16763079262169 16 293228729314067 15 

30240 - Kl * 2 + ^890~ KlK2 ~ 322560^ 2 + 90720 Kl + 60480 ^ 

75400838218273 , 4 - 599027907233393 . 14354418431047 , 2 , 

1890 + 3780 K ' Kl + 40 

544406144647673 , , , 49513057580797 , n , 6591403549631 « - 46629805507009 , , 

K\ X K% H K}°Kx + if' \Kl H K?K%+ 

1080 1 2 108 1 2 24 1 2 432 1 2 

58466928681791 7 9 730982464631 6 10 42641216914 5 n 973964951 4 12 

2160 ^ 2 + 180 KlK2 + 135 K ^ 2 + 120 K ^ 2 + 

139288171 3 451397 508229 15487 

T080— KlK ' -^890~ KlK2 + ^mim KlK2 " 72576^ 2 + 

5893118865913171 17 2424120194567633 16 29877830035592213 152 

13063680 1 + 181440 1 2 + 241920 1 2 + 

33355738800284053 14 , 85279687665510997 n 4 175231349967878723 12 5 

60480 + 60480 K ?«* + 77760 K "^ + 

91688471004674387 „ u , 21406776018978923 , n - 8051864567496943 „ „ 

K^KS H if^ifl H if? iff + 

38880 1 2 12960 1 2 10368 1 2 

12502971065814437 8 9 410022854780933 710 11705929429057 611 

51840 1 2 8640 1 2 2160 1 2 

7560435703027 t . 5t . 12 465465383273 , , 3 4996914619 . , 4 

iv n i\ 9 H i\i i\ 9 iVi iv 9 + 

25920 1 2 181440 1 2 36288 1 2 

1454610517 -2 ^ 15 415019,.. r _ 16 724301 , 7 17775777329026559 , 8 

if i if, ivi iv 9 H iv 9 H iv n + 

725760 1 2 22680 2 5225472 2 16329600 1 

59110336698456487 , 7 ^ 33916755588333461 , fi - 112823528586303187 r ^ 15 ~ 

ifi if 2 H K} b Kn H K} 5 KZ+ 

1632960 1 90720 1 2 60480 1 2 

121270297167625891 14 4 262383267925844563 13 5 446682230590568881 12 6 

22680 Xl ^ 2 + 27216 Kl K ' 2 + 38880 Kl * 2 + 

18047652078347753 n ^ 7 66470311691332679 r ^ 8 461244600701317 9 9 

1944 Al K2 + 12960 Kl K2 + 240 ^ 2 + 

2053612405901317 810 95019500270231 711 1801642236389 612 

4320 1 2 1296 12 288 12 

1068304780441 513 _ 356581981 414 _ 3715798969 315 27612313 216 

5040 1 2 216 1 2 60480 1 2 15120 1 2 

42818873 17 6111499 18 1697692411053976387 19 225292842432658127 18 

~ 1632960 1 2 + 32659200 2 + 653184000 1 + 2332800 1 2+ 

12073089792078645727 , 7 - 14877683079746729947 , fi - 35635183308116238361 , * , 

if i if o H K\ b Ki H if, 15 if 9 + 

10886400 1 2 2419200 1 2 1814400 1 2 

17990028013859879689 „ Ut ^ 24197213631925141831 ~ 6374234328455740363 , 2 - 

ifiifi 5 + kFkz + if, 12 ifZ+ 

453600 1 2 453600 1 2 129600 1 2 

906251536245123461 „ u „ 9035338464848376767 , n „ 85175951228570221 q , n 

if^iff H if, 10 if? H if? if 10 + 

28800 1 2 648000 1 2 20250 1 2 

23941624054521991 8 n 4310281756047181 7 12 5571935169950327 6 13 

28800 1 2 43200 1 2 907200 1 2 

23050830711047 . , 4 419715210991 , r ^ 15 48764703667 , , 6 

Ki Kn iv n ivn H iVi iv 9 + 

226800 1 2 129600 1 2 2419200 1 2 

4501185937^ 2 ^ 17 18779261^, 18 64742261 19 41816028466101527^ 20 

10886400 1 2 ~ 2041200 1 2 + 1306368000 2 + 6804000 1 + 

20744074405363158427 r ^ 19 ^ 3296051006124220457 r ^ 18 o 399315806473425521 , 7 , 

iVi iv 2 H iVi iv 9 H ill iv 9 + 

81648000 1 1020600 1 2 20160 1 2 

9084520466090723077 16 4 428717807069556189091 15 5 161475419175188528617 14 6 

1296W Al ^ 2 + 27216M K{ K ' 2 + 680400^ Ky * 2 + 
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593574184668411599 13 ^ 7 17599421384527454017 ^ 12 ^ 8 5041937738484379241 ^ n r ^ 9 
2400 Kl K2 + 97200 Kl * 2 + 54000 Kl * 2 

5438875343509708877 , , 1328632848231694729 Tjr% , x 20953128568172143 , , 2 
162000 1 2 162000 1 2 16200 1 2 

10792829938226569 7 13 1119984288812191 6 14 187714666081 s 15 
90720 1 2 + 226800 1 2 + 453600 1 2 + 

1090393669307 4 16 434088485513 317 2753717513 218 

453600 12 + 8164800 12 4082400 12 + 
562928201 iq 9330469 9n 206973837048951639371 21 

jy jy i y jy £ u _i_ t>^ z i _i_ 

81648000 1 2 163296000 2 " 14370048000 1 
53873357653341570569 20 40267035928596854137 19 2 4079484985047403119493 18 3 
81648000 1 2 + 4354560 1 2 + 65318400 1 2 ' 

379052188965132965269 17 4 13196982051156523866581 16 5 44245435623157935554441 15 6 
1555200 Xl ^ 2 + 21772800 Xl ^ 2 + 43545600 Kl * 2 

8616554448586239218387 t . 14t . 7 21437939380897068658541 r . 13r . 8 1133763862909311000887 r _ 12l ._ 9 
7257600 1 2 21772800 1 2 1944000 1 2 

39776679878739884231 n 10 15594884456933809459 10 n 15458199431760647 9 12 
162000 Kl ^ 2 + 216000 ^ ^ 2 + 1080 ^ 2 

173186300601405355 „ 890557841004984769 - , 4 42723412499535059 , i - 

if? if 13 h if; id 4 h k?ki 5 

96768 1 2 7257600 1 2 14515200 1 2 

-1259434901063897 516 17420157160631^ 4ifl7 774855209473 K 3 R is 

21772800 1 2 32659200 1 2 + 21772800 1 2 + 

-35002695431 , , 9 13528453 „ „ 20 108949189 , x 

ill iv 9 ' H iv i iv 9 iv 9 + ... 

46656000 1 2 1458000 2 1916006400 2 
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FIG. 1: The quantity S = Ki c (2; 0)/ifi c (2; R) — 1, which represents the reduced shift of the 
critical temperature of the system with anisotropy R from its 2D limit, is plotted vs P = R 2 / 3 . 
A continuous line interpolates our estimates, whose error bars are smaller than the width of the 
line, except for very small R. The dashed line (hardly visible except for small R ), which is 
superimposed to the continuous one, is the result of a fit of the expression f(R) = aR 9 + c to our 
data for 0.05 < R < 3.4. The values of the fit parameters are a « 1.245, g « 0.661 and c w 0.221. 



TABLE I: Estimates of the critical inverse-temperatures of the anisotropic system for various 
values of R. Only the two smallest-i? estimates are biased, the remaining ones being unbiased. 



R 


K lc (2;R) 


R 


K lc (2;R) 


R 


K lc (2;R) 


3.4 


0.13901(6) 


1.3 


0.20727(3) 


0.15 


0.3562(3) 


3.0 


0.14724(6) 


1.0 


0.22710(3) 


0.1 


0.3776(5) 


2.6 


0.15693(6) 


0.8 


0.24394(3) 


0.075 


0.391(2) 


2.2 


0.16861(5) 


0.6 


0.26537(3) 


0.05 


0.409(3) 


1.9 


0.17914(5) 


0.4 


0.29460(6) 


0.0125 


0.454(4) 


1.6 


0.19176(4) 


0.2 


0.3395(2) 


0.005 


0.473(4) 



TABLE II: Estimates, by first-order DAs, of the critical values G s of the normalized ratios 
G^iKi) as K x -»tfi c (2;0). 



critical value 


Gf 


gT 


Gf 


Gf 


of 


Gf 




1. 


1.00(1) 


1.01(1) 


1.00(2) 


1.02(8) 


0.95(15) 


critical value 


Gf 


Gf 


Gf 




Gf 






1.04(5) 


.99(4) 


0.99(1) 


0.99(1) 


.99(3) 


0.98(13) 
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FIG. 2: A blow-up of the lower left corner of FigfTJ showing the plot of the reduced critical- 
temperature shift S = Ki c (2; 0) / K\ c (2; R) — 1 vs P = i? 2 / 3 . The continuous line is a fit of the 
expression V/[ln(R/W)] 2 to our estimates of S (circles) for 0.005 < R < 0.15. The values of 
the parameters are V ~ 11.34 and W ~ 12.7. The short-dashed line is the result of a fit of the 
expression f(R) = aR 9 ' , with a « 1.08 and g' « 0.354, to the same set of data. The long-dashed 
line is the result of the same fit as in FigQ] of the expression f(R) = aR 9 + c, with parameters 
a 1.245, g rs 0.661 and c « 0.221, to all data for 0.05 < R < 3.4. The triangles are smalkR 
estimates of S taken from the simulation of Ref. 21] 



TABLE III: Estimates of the universal ratios Q s , computed by extrapolating first-order DAs of 
the effective ratios in eq. (j4*3]) . 



Qx 


Q2 


Q 3 


Qa 


Q 5 


Q& 


1.584(6) 


1.436(2) 


1.287(9) 


1.24(1) 


1.18(1) 


1.12(9) 
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0.0 0.10 0.20 0.30 



FIG. 3: A plot of SB = l/2K 2c (2; R) (circles) vs PB = i? 2 / 3 , in the small R region. A continuous 
line interpolates our estimates. 




FIG. 4: Estimates of the exponents j(2; R) (circles), v(2;R) (triangles) and of the ratios 
v(2; R)/ , y(2; R) (rhombs), plotted vs the anisotropy parameter R. All these quantities are com- 
puted by DAs biased with the critical temperature and they are normalized to our estimates of 
their values at R = 1. The horizontal lines are bands of 0.5% deviation from our estimates of 
7 (2; 1), !/(2;l) and 7 (2; l)/i/(2; 1). 
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0.0 0.2 0.4 0.6 0.8 1.0 

1 - KVK1c(2jR) 

FIG. 5: The effective exponent j e ff(2;Ki,R), computed by PAs, for the susceptibility of the 
anisotropic system is plotted vs r(2; R) = 1 — K\/K\ c {2\ i?) for various fixed values of R indicated 
on the corresponding curves. 




0.0 0.2 0.6 0.6 0.8 1.0 

1-K2/K2c(2;1ZR) 

FIG. 6: The effective exponent 7 e //(2; if 2, R), computed by PAs, for the susceptibility of the 
anisotropic system vs t(2;R) = 1 — K2/K2 C {2;R) for various fixed values of R indicated on the 
corresponding curves. 
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FIG. 7: The effective exponent j e ft(2; Ki, R) of the susceptibility, as obtained from the 
scaling form eq.(|l5"]). plotted vs t(2;R) = 1 — K\jK\ c {2\ R). From top to bottom R = 
0.0025,0.005,0.015,0.03,0.05,0.075,0.1,0.15. These values of R coincide with the eight smallest 
values used in Fig|5j 
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